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Abstract. We propose to explore the so-far poorly 
measured cosmic ray and gamma-ray sky (accelera- 
tor sky) in the energy range from 10 TeV to lEeV. 
New physics questions might be addressed in this 
last remaining observation window of gamma-ray 
astronomy. The very high beam-energies provided 
by Cosmic accelerators and the air-shower detection 
technique naturally imply an entanglement between 
fundamental questions of astroparticle physics and 
particle physics. The new large-area (10 km^) wide- 
angle (Isr) air Cherenkov detector SCORE (Study 
for a Cosmic ORigin Explorer) is based on non- 
imaging Cherenkov light-front sampling with sen- 
sitive large-area detector modules of the order of 
Im^. The lateral photon density and arrival- time 
distribution will be sampled up to large distances 
from the shower core. The physics motivations, the 
detector concept and first simulation results will be 
presented. 

Keywords: Non-imaging Cherenkov, Ultra-high En- 
ergy Gamma-Ray Astronomy, New Experiment 

I. The last Gamma-ray Observation Window 

The sensitivity level of existing and currently planned 
gamma-ray detectors is optimized for the very high- 
energy regime (VHE, 100 GeV < E < 30 TeV). The sen- 
sitivity to the ultra-high energy gamma-ray regime 
(UHE, E > 30 TeV) is limited because previously, the 
trend in development of detectors for gamma-ray as- 
tronomy was dominated by the focus on low energy 
thresholds. So far, the number of sources detected above 
10 TeV is less than 10 (see Figure[T]i, and above 100 TeV 
no detection has been reported so far A number of 
fundamental physics questions need to be addressed 
in the UHE gamma-ray regime. Future and planned 
experiments such as CTA fTH, HAWC lH?) or LHAASO 
Q will improve the situation in the UHE gamma-ray 
regime. However, due to dropping event statistics with 
rising energy, the key to UHE gamma-ray astronomy 
is a very large instrumented area of the order of 10 
to lOOkm^. While such large instrumented areas seem 
impracticable using the well-established imaging air 
Cherenkov technique (order of 10000 channels/km^), 
the non-imaging air Cherenkov technique provides a 
complementary possibility that comes along with some 
advantages, such as a small number of channels per km^ 
(few 100 per km^) and a wide field of view (order of 
sr). We have started the development of SCORE (Study 
for a Cosmic ORigin Explorer), a lOkm^ prototype 
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Fig. 1. The Ultra High-Energy Gamma-Ray Sky in Galactic coor- 
dinates. Shown ai'e sources so far detected with more than 5 a above 
E>10TeV. Status as of May 2009. 



detector for a ground-based wide-angle large-area air- 
shower detector for non-imaging gamma-ray astronomy 
and cosmic ray physics from 10 TeV to 1 EeV. With a 
wide field of view of the order of 1 sr (continuously 
monitoring a large part of the sky) and a focus on 
the highest energies, this project is complementary to 
(yet independent of) existing and planned experiments. 
SCORE is understood as a first step towards a lOOkm^ 
detector of similar design, H; SCORE (Hundred Square- 
km Cosmic Origin Explorer). 

In addition to its capabilities for gamma-ray observa- 
tions, a large area air-shower array such as SCORE is the 
ideal detector for cosmic ray physics, allowing observa- 
tions of primary cosmic rays in the domain of transition 
between the Galactic and the extra-galactic component. 
The detection technique of SCORE is based on shower 
front sampling and timing with Cherenkov light (non- 
imaging technique). SCORE provides a unique approach 
to measure the entire longitudinal development using the 
shower-front arrival-time distribution (at distances from 
the shower core >100m), allowing detailed spectral 
and composition measurements, and allowing gamma- 
hadron separation via reconstruction of the shower- 
depth. With an instrumented area of 10 km^, SCORE 
covers an energy range for measurements of the primary 
cosmic ray spectrum and the chemical composition 
from 100 TeV - overlapping with direct cosmic ray 
observations - to 1 EeV0- overlapping with fluorescence 
measurements. 

In the next section the main physics objectives of 

' The number of detected cosmic rays above 1 EeV is expected to 
be 10 per year 
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SCORE will be addressed, followed by a description of 
the detector principle ( section HiHi and first simulation re- 
sults (section ITVll. The event reconstruction is described 
in a dedicated contribution to this conference ifTTI . 

II. Physics Objectives 
A. Cosmic Ray Origin 

The energy range from 10 TeV to several 100 TeV 
is the most natural energy regime for a solution of the 
almost 100-year-old mystery of the origin of cosmic 
rays. Gamma-rays are emitted at the site of cosmic 
ray acceleration via hadronic interactions of the cos- 
mic rays with ambient matter, producing neutral pions 
subsequently decaying into gamma-rays (and neutrinos 
from charged pions). At least some Galactic accelerators 
- the pevatrons - are expected to accelerate cosmic 
rays via shockfront acceleration [91 up to knee energies 
(Ri4PeV). The gamma-ray spectral energy distribution 
of a pevatron is expected to reach up to several lOOTeVH 
(see Figure O. An additional spectral hardening (flatten- 
ing) below the cut-off energy is predicted in diffusive 
shock acceleration theory when non-linear effects are 
taken into account llT9l , 

In the VHE regime, an ambiguity in gamma-ray 
spectra arises from the fact that gamma-rays can also 
be produced by leptonic acceleration processes (inverse 
Compton scattering). Therefore, if observations only 
provide spectra up to few tens of TeV, the underlying 
production mechanism of a gamma-ray signal needs to 
be resolved with precise spectral measurements over a 
wide range in energy (multi-wavelength) and detailed 
comparisons to predictions from hadronic and leptonic 
scenarios. Luckily, the ambiguity of gamma-ray spectra 
disappears in the UHE regime, where the inverse Comp- 
ton scattering cross-section enters the Klein-Nishina 
regime, loosing efficiency. This results in a softening (i.e. 
steepening) of gamma-ray spectra from leptonic acceler- 
ators beyond 10 TeV. As opposed to that, a hard (i.e. flat) 
gamma-ray spectrum continuing up to few hundred TeV 
would be a clear signature of hadronic acceleration. For 
comparison, the observed spectral energy distribution of 
the shell-type supernova remnant RXJ 1713.7-3946 |1 1 is 
shown as a dotted line in Figure |3] The spectrum of this 
supernova remnant - as well as of any other object so far 
detected in the VHE regime - can still be accomodated 
within leptonic acceleration scenarios. Comparing the 
mesured spectra with expected pevatron spectra shows 
that an expectedly existing population of PeV acceler- 
ators might have been missed out completely by VHE 
experiments, even more so when considering that such 
accelerators might be widely extended structures larger 
than the field of view of the most sensitive instruments 
(lACTs, 3-5°). Recent detections of partly unidentified 
sources (H.E.S.S., Milagro) yield hard gamma-ray spec- 
tra that extend up to few tens of TeV. These results 

^The cut-off energy in gamma-ray spectra is roughly a factor of 10 
lower than in the underlying hadronic spectrum. 



underline the necessity to carry out observations in 
the UHE regime. As opposed to previous experiments, 
SCORE is specifically designed for the UHE gamma- 
ray regime and for spatially extended emission. With the 
SCORE experiment, we aim at identifying cosmic ray 
accelerators by searching for hard spectra reaching up to 
few 100 TeV, therwith contributing to the resolution of 
the mystery of the origin of cosmic rays. The observation 
of UHE gamma-ray emission from extended gas clouds 
illuminated by nearby cosmic ray sources, might open up 
a possibility to map the Galactic cosmic ray acceleration 
efficiency. 

B. Propagation and Absorption 

At such high energies as considered here, gamma-rays 
are attenuated via e+e^ pair production with the photons 
of the Galactic interstellar radiation fields (IRF) and 
with the Cosmic Microwave Background (CMB). It was 
shown that for Galactic objects the attenuation reaches a 
maximum around 100 TeV from the Galactic IRF and at 
2 PeV from the CMB |16|. While the former depends 
on the Galactic longitude and local radiation fields, 
the latter is universal. If the distance of the observed 
gamma-ray sources is known, the density of the IRF 
might be inferred from the strength of the absorption 
by pair production (or from spectral features). Inversely, 
Galactic absorption might also open up a new possibility 
to infer distances from the measurement of gamma-ray 
spectra, if the IRF in the line of sight is known. Such a 
new method for distance estimation of Galactic objects 
might also be possible if the universal absorption by 
the CMB could be measured. The expected attenuation 
by pair production might be altered by photon/axion 
conversion (e.g. |18|). Photons produced at the source 
travelling through the Galactic magnetic field might 
convert into axions propagating without absorption. If a 
reconversion of these axions back into photons happens 
before arriving at earth, the photon signal would appear 
to be stronger than expected. The same effect could arise 
if photon/hidden-photon oscillations ll20l would occur 
Another effect that might alter the expected absorption 
by pair production is the modification of the e+e^ 
pair production threshold in case of Lorentz invariance 
violation. 

C. The Local Supercluster 

Ultra-high-energy cosmic rays from the local super- 
cluster are expected to interact with the CMB, initiating 
intergalactic pair cascades that could be observable in 
gamma-rays. The emission is expected to match the 
large scale of the local super-cluster structure and would 
be measurable as anisotropic diffuse emission in the 
total field of view of the experiment. Alternativeley, the 
accelerators of ultra-high-energy cosmic rays might ex- 
hibit point-like emission or halo-like structures resulting 
from the interactions of the accelerated particles with 
the surroundings of the source. 
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D. Cosmic Ray Composition 

By measuring the longitudinal development of 
hadronic air-showers (using the shower-front arrival-time 
distribution), detailed spectral and cosmic ray compo- 
sition measurements are possible from 100 TeV to 1 
EeV. SCORE will detect 2 x 10^ cosmic rays per year 
above 1 PeV and 10 cosmic rays per year above 1 EeV. 
Furthermore, the reconstruction of the shower depth 
will allow an independent measurement of the inelastic 
proton-proton cross-section overlapping and exceeding 
LHC energies. 

III. The Detector Concept 

The detector principle is based on the shower-front 
sampling technique using Cherenkov light. A large array 
of wide-angle light-sensitive detector stations is used. 
The concept of the detector modules used as working as- 
sumption for the SCORE detector was adapted from pre- 
vious gamma-ray experiments, such as THEMISTOCLE 
ifTOll . 13 1, HEGRA AIROBICC |15|, or BLANC A |8|. 
Similar detector modules are also used in the currently 
being built TUNKA array for cosmic-ray physics ||6l. As 
compared to the previous experiments, three aspects of 
SCORE will be different: an instrumented area larger by 
more than an order of magnitude, larger detector station 
areas and larger inter-station spacing (see Table |l]i. Ad- 
vances made in technology allow improvements to the 
original detector components such as improved photo- 
sensitive detectors, fast trigger and readout electronics, 
therewith allowing a measurement of the Cherenkov 
photon arrival time distribution. 

A very large instrumented area (lOkm'^) is required 
to reach sufficiently large event statistics at ultra-high 
energies, and is achieved with a low array density, i.e. 
large inter-station spacings. A reasonable value for the 
detector station spacing can be derived from the lateral 
photon density function (LDF) of Cherenkov light at 
observation level, shown in Figure |2] Within a radius 
of 120 m around the shower core position the LDF 
is roughly constant, but shows a large spread from 
shower to shower. Fluctuations are much lower beyond 
120 m. With the envisaged station spacing of 1 00-200 m, 
SCORE will primarily be sensitive to the outer part 
of the LDF (promising good reconstruction quality), as 
illustrated with the solid lines in Figure |2l The low 
photon density far away from the shower core justifies 
the chosen large individual detector station areas. For 
comparison, the corresponding sensitive range of the 
AIROBICC experiment is also shown (dashed lines), 
illustrating the novel approach chosen with SCORE. 
Another important aspect for the reconstruction with 
the SCORE detector is the usage of the full timing 
information from the arrival time distribution of the 
Cherenkov photons at each detector station. The event 
reconstruction of SCORE is based on the combination of 
information from the lateral photon density distribution 
and the arrival time distribution of Cherenkov photons 
(see ifTTl . this conference). 
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Fig. 2. Lateral photon density function (LDF) of Clierenkov light 
on sea level for airshowers initiated by a 10 TeV gamma-ray shower 
(grey area) and a 100 TeV gamma-ray shower. The sensitivity level 
of one SCORE detector station is indicated by the soHd line. For 
comparison, the corresponding sensitivity level for AIROBICC is also 
shown (dashed grey box). 

Each detector station consists of four photomulti- 
plier tubes (PMTs) equipped with four light-collecting 
Winston cones of 30° half-opening angle pointing to 
the zenith. All four modules (PMT-i-cone), including 
the trigger, readout electronics and communication are 
planned to be encased in a box equipped with a sliding 
lid. The advantages of using four PMT channels per 
station are the possibility to suppress false triggers from 
nightsky background (NSB) light by a local coincidence 
trigger condition and the requirement of a large light 
collecting area (a). A total area of a = 0.5 m^ can be 
achieved when using four 8" PMTs and a Winston cone 
height of 0.5 m. The final design of the Winston cones 
is work in progress. Different solutions are envisaged 
for the cone structure (plastic mold, segmented metal 
sheets) and for the reflective inner surface (vaporization, 
reflective foils). A fast signal readout and digitization in 
the GHz regime are needed. Different solutions such as 
analog ring samplers or domino ring samplers (DRS) are 
under study. We are currently testingusage of the DRS 4 
chip that was developed by the PStI and which is also 
used by the MAGIC experiment. 

The envisaged trigger scheme encompasses two lev- 
els. A local station trigger is issued when the sum of 
all four PMT signals passes a given threshold (few a 
above NSB level). In the second level, a two-fold station 
coincidence is required from a sub-array of detector 
stations (3 x 3 or 4x4) and the signals are then sent to the 
data acquisition. Individual sub-arrays overlap to avoid 
dead trigger zones. Alternatively a smart station concept 
could be implemented, issuing a trigger based on a next- 
neighbour station condition. For a large array such as 
SCORE, with large inter-station spacings, relative timing 
is a critical issue. A relative timing accuracy of better 
than 5 ns is required to limit the impact on reconstruction 
(mainly angular resolution). For time synchronization we 
plan to use optical fibers for clock distribution to the 

^http://midas.psi.ch/drs| 
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TABLE I 

Basic design characteristics of the SCORE detector and for the possible extension HiSCORE. Corresponding values 

FOR TUNKA, AIROBICC AND BLANCA ARE ALSO GIVEN. 





SCORE 


Hi SCORE 


TUNKA 


BLANCA 


AIROBICC 


THEMISTOCLE 


Instrumented area Affcm''] 


10 


100 


1 


0.2 


0.04 


0.08 


Detector station area a[m^] 


0.5-1.5 


0.5-1.5 


0.2 


0.1 


0.13 


0.5 


Field of View FoV[sterad] 


0.84 


0.84 


1.8 


0.12 


1 


0.008 


Inter- station distance d[m] 


100-200 


100-200 


85 


35 


30 


50-100 


Number of detector stations 


O(IOO) 


2601 


133 


144 


49 


18 




energy [TeV] 

Fig. 3. Sensitivity of SCORE (solid line) after 5 years of operation. 
Further sensitivities are given for H.E.S.S., CTA and HAWC. For com- 
parison, the expected pevatron spectrum (exponential cut-off energy 
was 4PeV, linear acceleration: solid grey line, non-linear acceleration: 
dashed grey line) and a fit to H.E.S.S. data from the shell-type SNR 
RXJ 1713.7-3946 (dotted grey line), as well as an upper limit by the 
KASCADE experiment 1 2 | are given. 



100 TeV and significantly extends the energy range in 
the UHE gamma-ray regime. 

V. Conclusion and Outlook 

With the developement of the SCORE detector we aim 
at opening up the last remaining observation window 
of gamma-ray astronomy, the UHE gamma-ray regime. 
Furthermore, cosmic ray spectral and chemical com- 
position measurements from 100 TeV to 1 EeV will be 
possible with SCORE. First simulation results show a 
competitive sensitivity above 100 TeV and demonstrate 
a good potential for new physics in the UHE regime. 

With further improvements of reconstruction and 
gamma/hadron separation power, we hope to improve 
the expected sensitivity of the SCORE detector With a 
larger array of lOOkm^ (HiSCORE, thick dashed line 
Figure O, the sensitivity could further be enhanced. 



local station clocks. 

An extension of the concept we studied so far is to 
equip the underside of each sliding lid of the station 
boxes with scintillator material. With such a setup, 
continuation of data taking during daytime would be 
possible using SCORE as a very large charged particle 
shower front sampling array. 

IV. First Simulation Results 

Air showers were simulated using CORSIKAv675 
lfT2l . The detector simulation {sim_score) was imple- 
mented on the basis the iact package |[5|. We included 
signal pulse shaping and used basic assumptions on PMT 
quantum efficiency and Winston cone reflectivity. In or- 
der to estimate the sensitivity of the SCORE detector, the 
effective trigger area was calculated for mono-energetic 
gamma-rays and protons. These results were used to 
estimate the number of signal events (for a pevatron- 
like spectrum as in Figure [3]) and background events 
based on the parametrization of the cosmic ray spectrum 
by ifTSl . An energy-dependent angular resolution as 
found in simulations (see [11], this conference) and a 
Quality factor for gamma-hadron separation of 1.5 were 
assumed. The point-source sensitivity of SCORE was 
calculated for a detection of at least 50 gamma-rays 
at a significance level of 5 cr in 5 years (see Figure |3] 
thick solid line). The SCORE sensitivity is competitive 
with existing and currently planned instruments around 
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